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ABSTRACT 

We derive a simple, fast one-dimensional integral for the equatorial rotation curve of a thin disk with 
surface density modeled as a spheroid with axis ratio q. The result is simpler than standard 

expressions even in the limit of an infinitely thin disk (g ^ 0). 

Subject headings: galaxies: kinematics and dynamics 



L INTRODUCTION 

Existing expressions for the calculation of the equato- 
rial rotation curves of thin disks have undesirable fea- 
tures. Derivations using elliptic integrals are complicated 
and contain a singularity in the plane which is usually 
removed by evaluating the rotation curve slightly above 
the plane (e.g. Binney & Trcmaine 1984). Derivations 
in terms of Hankel (Toomre 1962) or Stieltjes (Evans & 
de Zeeuw 1992) transforms are poorly suited for numeri- 
cal calculation. Standard expressions for finite thickness 
disks require intrinsically two-dimensional numerical in- 
tegrals (e.g. Casertano 1983, Cuddeford 1993) or multi- 
dimensional tables (e.g. Dehnen & Binney 19981 ). 

In conducting some experiments on the compatibility 
of central dark matter cusps with galaxy rotation curves 
(see, e.g., Moore 2001), we needed a faster means of com- 
puting equatorial rotation curves including the effects of 
a finite disk thickness. By faster we mean a simple, non- 
singular, non-oscillatory, one-dimensional numerical inte- 
gral for the rotation curve of a disk with an arbitrary 
surface density profile S(i?). The solution was found 
by considering the equatorial rotation curves of flattened 
spheroids. Modeling disks as flat spheroids is little used, 
although the mathematical approach was developed by 
Brandt ( |l960D , Brandt fc Belt on ( |l962| ), Mes tel (|1963D , 
Lynden-Bel l fc P ineault ( |1978| ), Lynden-BeU ( |l989| ) and 
Cuddeford (|l993| ). These studies show that the equatorial 
rotation curve can be expressed as an integral over the 
Abel transform of the disk surface density, as also outlined 
in Binney & Tremaine (1984). In this paper we follow the 
same development, but then show that the resulting two- 
dimensional integral can be reduced to a one-dimensional 
integral with convenient properties for numerical calcula- 
tion. We will refer to our solution as the rotation curve of a 
spheroidal disk, since it is the rotation curve of a spheroid 
{p{B? + /q^)) defined by the surface density S(-R) of a 
disk. 

2. THE EQUATORIAL ROTATION CURVE OF A 
SPHEROIDAL DISK 

We develop our result follow ing the procedur es out lined 
in Binney & Tremaine (|l984|) and Cuddeford ( |l993|) , al- 
though we focus only on the equatorial rotation curve 



rather than the global potential. For a density distribu- 
tion p{m?) with m? ~ + /q^ , the equatorial rotation 
curve is simply 
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where the density distribution is the Abel transform of the 
surface density 
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for a sufficiently regular surface density. For example, the 
Abel transform of an exponential disk, E = Sq exp(— i?/_Rc;), 
is 
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where Kq(x) is a modified Bessel function, and the rotation 
curve is 
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Cuddeford (|T99|). For g = we obtain the standard ana- 
lytic result by Freeman ( |197Q| ). 

Our extension to these results is that we can reduce the 
calculation to a one-dimensional integral by reversing the 
order of integration to find that 

vl{R) = du— [i?(<^, a) - F{<l>, a)] (5) 

where = (1 — q^)v? / R^ and sin0 = min(l, i?/u)|^ The 
integrand has an integrable, logarithmic singularity at a = 
1 {u = R) for an infinitely thin disk (5 ^ 0). The result 
simplifies further to 
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^The special functions FU,a) = ff de{l - a^sin^e)-!/^ and 

E(ct>,a) = rfe{l - Q^sin^e)^/^ are elliptic integrals of the first 
tf 

and second kind respectively. 



if we use the Carlson forms for elliptic integrals^ rather 
than the Legendre forms. The function R]j{x,y, z) has a 
logarithmic singularity when x = y = 0, again correspond- 
ing to a = 1 {u = R) with q = 0. Integrating by parts, to 
convert dT,/du into E, and taking the limit g — > leads to 
the standard expression for the circular velocity in terms 
of elliptic integrals evaluated in the disk plane (e.g. Bin- 



ney & Tremaine 1984). Integrands based on (E/dR, like 
eqns. (5) and (6), have a weaker singularity than those 
based on and the singularity is easily removed by 

giving the disk a non-zero axis ratio q. 

3. SUMMARY 

Our expressions for the equatorial rotation curves of 
spheroidal disks of surface density and axis ratio q 

are simple (even compared to standard expressions for in- 
finitely thin disks) and can be computed very rapidly due 
to the excellent numerical properties of the Carlson ellip- 
tic integrals. The one disadvantage of the result is that 



real disks (e.g. de Grijs (199?) and references therein) 



are better modeled as separable functions of R and z (e.g. 
E(i?) exp(-|z|/_ff)) than as spheroids {p{R^ + z"^ / q^)) with 
axis ratio q. This distinction will be important in some cir- 
cumstances (e.g. detailed studies of vertical equilibrium in 
stellar disks), but should matter little for the global rota- 
tion curve if the axis ratio is adjusted to roughly match 
the central disk scale height. In any case, the simplicity of 
our final expression makes it useful even for infinitely thin 
disks. 
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■^Carlson's elliptic integral of the second kind is defined by 

-RD(a;,2/,2) = l/p (;^ + ^)l/2(;^+t,'^)l/2(;^ + ^)3/2 ' 



